Thermodynamics of quantum degenerate gases in optical lattices 
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The entropy-temperature curves are calculated for non-interacting Bose and Fermi gases in a 3D optical 
lattice. These curves facilitate understanding of how adiabatic changes in the lattice depth affect the temperature, 
and we demonstrate regimes where the atomic sample can be significantly heated or cooled by the loading 
process. We assess the effects of interactions on a Bose gas in a deep optical lattice, and show that interactions 
ultimately limit the extent of cooling that can occur during lattice loading. 
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I. INTRODUCTION 

Tremendous advances have been made in the preparation 
and control of bosonic and fermionic atoms in optical lattices 
(e.g. see QESSHSSIES OH)- In experiments the 
gas is typically prepared in the optical lattice by a slow load- 
ing procedure that begins with a weakly trapped gas and no 
lattice. During loading, the lattice is turned on in some pre- 
scribed way, and the atoms are localized into the tightly con- 
fining potential wells of the optical lattice. This process is 
accompanied by a massive redistribution of the energy states 
of the system, and it is poorly understood how the loading 
process affects the properties of the atoms, such as their tem- 
perature. Many of the physical phenomenon that are suitable 
to experimental investigation in optical lattices are sensitive to 
temperature and it is therefore of great interest to understand 





Figure 1: Schematic diagram of system under consideration: a col- 
lection of atoms in a uniform lattice of controllable depth. 



how the temperature of a quantum degenerate gas changes 
with lattice depth. Experimental results by Kastberg et al. 
fllll in 1995 showed that loading laser cooled atoms into a 
three-dimensional optical lattice caused the atoms to increase 
their temperature 12311 . Recent studies have shown that there 
is rich range of behaviour that can be expected to occur dur- 
ing the loading process at temperatures much lower than those 
explored by Kastberg (e.g. see OfijlfT^ll V 

In this work we compare and contrast the behaviour of 
Bose and Fermi systems in optical lattices as a function of 
the lattice depth. While some work has been undertaken in 
the case where an external harmonic potential is also present 
(but restricted to a non-interacting tight-binding approxima- 
tion) n, here we will restrict our attention to the uniform lat- 
tice. Schematically, our system is shown in Fig. ^ a system 
of atoms confined in a 3D cubic lattice. As the lattice depth 
increases the potential changes from that of a uniform box 
potential to that of a deep lattice potential. The fundamental 
question we wish to address is how the properties of the equi- 
librium state, in particular the temperature, change under this 
loading procedure. To identify the final thermodynamic state 
of the system we assume that the loading procedure is isen- 
tropic. In practice experiments appear to be approximately 
reversible when the loading is performed slowly, as has been 
investigated in Ref. rial . 

The paper is organized as follows: In section [I]] we intro- 
duce the theoretical approach we use for calculating the single 
particle spectrum and thermodynamic properties of an ideal 
quantum gas in an optical lattice. The results of this formalism 
are presented and discussed in section [TO] In section HV1 we 
address the effects of interactions in application to a bosonic 
system, before concluding. 



II. IDEAL GAS FORMALISM 

The single particle spectrum completely determines the 
thermodynamic properties of an ideal gas. We consider a 
cubic 3D optical lattice made from 3 independent sets of 
counter-propagating laser fields of wavelength A, giving rise 
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to a potential of the form 

V 

VL a tt(r) = — [cos(2fcr) + cos(2ky) + cos(2kz)}, (1) 

where k = 2ir/\ is the single photon wavevector, and V is the 
lattice depth. We take the lattice to be of finite extent with a 
total of N s sites, consisting of an equal number of sites along 
each of the spatial directions with periodic boundary condi- 
tions. The single particle energies e q are determined by solv- 
ing the Schrodinger equation 
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(2) 



for the Bloch states, i/)q(r), of the lattice. For notational sim- 
plicity we choose to work in the extended zone scheme where 
q specifies both the quasimomentum and band index of the 
state under consideration [24]. By using the single photon 
recoil energy, Er = h 2 k 2 /2m, as our unit of energy, the en- 
ergy states of the system are completely specified by the lattice 
depth V and the number of lattice sites N s (i.e. in recoil units 
e q is independent of k). 

For completeness we briefly review some important fea- 
tures of the band structure of Eq. (|2) relevant to the ther- 
modynamic properties of the system. The smoothed density 
of states for the system for various lattice depths is shown in 
Figs. |2ja)-(d). For sufficiently deep lattices an energy gap, 
e gap , will separate the ground and first excited bands (see Fig. 
0c)). For the cubic lattice we consider here, a finite gap ap- 
pears at a lattice depth of V « 2Er [25] (marked by the ver- 
tical asymptote of the dashed line in Fig. |2je)). For lattice 
depths greater than this, the gap increases with lattice depth. 
In forming the gap, higher energy bands are shifted upwards in 
energy, and the ground band becomes compressed — a feature 
characteristic of the reduced tunneling between lattice sites. 
We refer to the energy range over which the ground band ex- 
tends as the (ground) band width, cbw (see Fig. EJc)). As is 
apparent in Fig. EJe), the ground band width decreases ex- 
ponentially with V, causing the ground band to have an ex- 
tremely high density of states for deep lattices. 

Our primary interest lies in understanding the process of 
adiabatically loading a system of N p bosons or fermions into 
a lattice. Under the assumption of adiabaticity the entropy re- 
mains constant throughout this process and the most useful 
information can be obtained from knowing how the entropy 
depends on the other parameters of the system. In the thermo- 
dynamic limit, where N s — > oc and N p — ► oo while the filling 
factor n = N p /N s remains constant, the entropy per particle 
is completely specified by the intensive parameters T, V, and 
n. The calculations we present in this paper are for finite size 
systems, that are sufficiently large to approximate the thermo- 
dynamic limit. We would like to emphasize the remarkable 
feature of optical lattices that V is an adjustable parameter, in 
contrast to solid state systems where the lattice parameters are 
determined by the constituent atoms and are immutable. 

We determine the entropy as follows: We calculate the sin- 
gle particle spectrum {e q } for given values of N s and V. 
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Figure 2: (a)-(d) Density of states for a N s « 3 x 10 4 site cubic 
lattice at various depths. For a depth of approximately V ~ 2Er a 
gap develops in the density of states. In (c) we illustrate the energy 
gap e gap and ground band width taw- Points are determined by 
numerically averaging the exact spectrum over a small energy range, 
(e) The dependence of the energy gap (e gap dashed line) and ground 
band width (ebw solid line) on the lattice depth (see the text). 



Working in the grand canonical ensemble we evaluate the 
grand partition function Z, according to 



log2 = ±J>g(l±. 



-/3(e q -M) 



(3) 



where /i is found by ensuring particle conservation, and F 
(B) refer to the case of fermions (bosons). The entropy of the 
system can then be expressed as 



S = k B {\ogZ + (3E~ix(3N p ). 



(4) 



where (3 = 1/kgT, and E = — d In Z /d(3 is the mean energy. 

Most current Fermi gas experiments are realized with a 
mixture of two internal states. This is required because s- 
wave elastic collisions, needed for re-equilibration, are pro- 
hibited by the Pauli principle for spin polarized samples 12611 . 
The theory we present here is for the spin polarized case, but is 
trivially extensible to multiple components if the lattice poten- 
tial is spin independent and the number of atoms in each com- 
ponent is the same: in this case all extensive parameters are 
doubled (e.g. {E, S}) and intensive parameters (e.g. {T, /i}) 
remain the same. The inclusion of interaction effects, which 
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will be important in the multiple component case, is beyond 
the scope of this paper. 



III. IDEAL GAS RESULTS 
A. Effect of lattice loading on temperature 



the degree of condensation will be greatly affected by adia- 
batic lattice loading. For instance, consider the adiabatic pro- 
cess indicated by the dashed line and labeled C in Fig. 0b). 
The system starts as a Bose-condensed gas of free particles. 
However, as the lattice depth increases the condensate frac- 
tion decreases until the system passes through the transition 
point and becomes uncondensed. 



In Fig. [5]we show entropy-temperature curves for various 
lattice depths and filling factors n. These curves have been 
calculated for a lattice with 31 lattice sites along each spatial 
dimension, i.e. N s » 3 x 10 4 . A general feature of these 
curves is the distinct separation of regions where adiabatic 
loading causes the temperature of the sample to increase or 
decrease, which we will refer to as the regions of heating and 
cooling respectively (e.g. see Figs. |3](a)-(c) and (e)). These 
regions are separated by a value of entropy, 5*0, at which the 
curves plateau, and we note that this feature is more prominent 
on the curves for larger lattice depths. This plateau entropy is 
indicated by a horizontal dashed line and is discussed below. 
For the case of fermions with unit filling, shown in Fig. |3jf), 
this plateau occurs at So = 0, and only a heating region is 
observed. 

We now explicitly demonstrate the temperature changes 
that occur during adiabatic loading using two possible adia- 
batic processes labeled A and B, and marked as dotted lines 
in Fig. |3jc). Process A begins with a gas of free particles in 
a state with an entropy value lying above the plateau entropy. 
As the gas is loaded into the lattice the process line indicates 
that the temperature increases rapidly with the lattice depth. 
Conversely process B begins with a gas of free particles in a 
state with entropy below the plateau. For this case adiabatic 
lattice loading causes a rapid decrease in temperature. This 
behavior can be qualitatively understood in terms of the mod- 
ifications the lattice makes to the energy states of the system. 
As is apparent in Fig. EJe), the ground band rapidly flattens 
for increasing lattice depth causing the density of states to be 
more densely compressed at lower energies. Thus in the lat- 
tice all these states can be occupied at a much lower temper- 
ature than for the free particle case. As we discuss below, 
for both Bose and Fermi systems, Sq is the maximum entropy 
available from only accessing states of the lowest band. If 
S < So, the temperature of the system must decrease with 
increasing lattice depth to remain at constant entropy. Alter- 
natively, for S > So the occupation of states in higher bands 
is important, and as the lattice depth and hence e gap increases, 
the temperature must increase for these excited states to re- 
main accessible. 



1. Bosonic systems 

The temperature and entropy at which bosons condense 
generally changes with lattice depth, and is indicated by cir- 
cles on the S-T curves in Figs. [3Ja) and (b). We note that for 
high filling factors the condensation points for different lat- 
tice depths occur over a wide entropy range, suggesting that 



2. Fermionic systems 

In addition to the effect that lattice loading has on the abso- 
lute temperature of a Fermi-gas, it is of considerable interest 
to understand how the ratio of temperature to the Fermi tem- 
perature (Tp) l27il changes. Indeed, the ratio T/Tp is the 
standard figure of merit used to quantify the degeneracy of di- 
lute Fermi gases. In Fig. |5Jd) we show how T/Tp changes 
with adiabatic lattice loading for the same parameters used in 
Fig- EJc). This result indicates the typical behaviour seen: 
Below the entropy plateau where cooling is observed (e.g. see 
Figs. |3jc) and (e)), the ratio of T/Tp remains approximately 
constant, so that there is little change in the degeneracy of the 
gas. Above the entropy plateau where heating was observed, 
the ratio of T /Tp rapidly increases, so that in this regime the 
gas will rapidly become non-degenerate as it is loaded into 
the lattice. For the unit filled Fermi case (Fig. |3Jf)X there 
is no cooling regime, and heating is accompanied by a rapid 
increase in T/Tp for all initial conditions of the gas. 

We note that for n > 1 it is possible to observe a reduction 
in T/Tp for Fermi systems during the loading process. This 
occurs because Tp increases, because the Fermi energy lies in 
the excited band. We do not consider this case here and refer 
to the reader to Ref.l 131 for details. 



B. Entropy plateau 

In many S-T curves a plateau in the entropy is apparent. 
This occurs when the gap in the energy spectrum between the 
states of the ground and first excited bands is large compared 
to the energy width of the ground band, i.e. e gap > eew- 
In this case, there is an intermediate temperature range, suffi- 
ciently hot that all the states in the ground band are accessed, 
yet not hot enough for states in the next band to be accessed. 
Within the temperature range satisfying these conditions, the 
entropy remains approximately constant at the value corre- 
sponding to the saturated ground band contribution to the en- 
tropy - we refer to this value as the plateau entropy. We now 
provide analytic expressions for this plateau. 

The total number of microstates accessible to the ground 
band, f2o, for the cases of fermionic and bosonic particles is 
given by 
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Figure 3: Entropy versus temperature curves for a N s m 3 x 10 4 site cubic lattice, at various depths V = to 20Er (with a spacing of 2Er 
between each curve). Cases considered are (a) bosons with n = 1.0, (b) bosons with n = 4.0, (c) fermions with n = 0.25, (d) fermions 
(reduced temperature) with n = 0.25 (e) fermions with n = 0.8, and (f) fermions with n = 1.0. The entropy plateau is shown as a dashed line. 
The processes indicated by the paths labeled A,B and C are discussed in the text. For the bosonic systems the critical point for condensation 
on each curve is indicated with a hollow circle. 
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respectively! 2§| . The corresponding value of entropy S 
k B log Qq, i.e. the plateau entropy, is given by 



N p k B 



--log(l-n)+log(--l 



log 1 



1 



1 



log (1 + n) 



, (7) 
(8) 



which we have explicitly written in terms of the filling fac- 
tor, n = N p /N s , with the additional validity conditions 
1 < N p < N s (Fermions) and N Sl N p > 1 (Bosons). We 
note an important case for which the above approximation is 
invalid is for the Fermi system with Np = N s , i.e. when we 
have a filling factor of n — 1, for which So = 0. This case 
corresponds to the unit filling factor result shown in Fig. |3jf) 
where, as a result of the entropy plateau occurring at So = 0, 
only a heating region is observed. 
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Figure 4: Tf versus T for final lattice depths of (dashed line) V = 
10E R , (dash-dot line) V = 20E R , and (solid line) V = 30E R . 
These results are for the case n = 1 and for the parameters of 87 Rb. 
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IV. EFFECTS OF INTERACTIONS IN THE BOSONIC 
SYSTEM 

We now briefly comment on the role of interactions on the 
properties of the Bose gas loaded into an optical lattice. For 
situations where the number of bosons is commensurate with 
the number of lattice sites, as the lattice depth increases the 
system will eventually enter the Mott-insulating state |y, Hal . 
In this state the system exhibits a gapped excitation spectrum 
which is poorly described by the (gapless) non-interacting 
spectrum. Here we develop an analytic description for the 
unit filled system in the strongly interacting regime to assess 
the effects that interactions have on adiabatic loading, and the 
behavior of temperature in the deep lattice limit. We also refer 
the reader to related results in Ref . 11711 . 

In a lattice of depth greater than a few recoils, the system is 
well-described by the Bose-Hubbard model flal 



H 



BH 



<ij> J 



(9) 



where aj is the bosonic annihilation operator of a particle at 
site j = {j x , jy, jz}, and the sum is over nearest neigh- 
bouring lattice sites. The interaction parameter U and the tun- 
neling parameter J can be determined by band structure cal- 
culations fl8il . 

In a deep lattice the tunneling parameter is exponentially 
suppressed, and can be taken to be approximately zero (2^1 . 
In this limit Fock states of the aj operators diagonalize the 
Bose-Hubbard Hamiltonian, and in the low temperature limit 
k B T < U and for < n < 1 we find that 



/' 



U 



(n-l+ y/3n{2 - n) + lj 



2(2 - n) 



(10) 



For the unit filled lattice, n = 1, the chemical potential is 
approximately /i = U/2 and one can derive an analytic ex- 



pression for the grand canonical partition function given by 
(see Ref. Q for details) 



(l + tf 3 (0,e-^) 



(ID 



where ^3(2,(7) = 1 = 2^^L 1 <7" cos(2nz) is the Elliptic 
Theta Function. From these results we can determine the en- 
tropy of the system using Eq. @. 

To assess the effect that interactions have on loading, we 
use entropy comparison between the initial system of zero 
lattice depth (i.e. the uniform system) and the final deep 
lattice case described by Eqs. dlOt and il Q . For the uni- 
form case we describe the system using the Bogoliubov ap- 
proximation O, in which case the quasiparticle spectrum 

for the uniform system is given by e q = )J(£q) 2 + 2itpe^, 
where = fi, 2 q 2 /2m, p — n/(A/2) 3 is the density and 
u = 4ira s h 2 /m, with a s the s-wave scattering length. As- 
suming that the quasiparticle occupation is given by the Bose 
distribution /(e q ) = [e^ £q — the entropy can be numer- 
ically evaluated according to 



Su„i = k B £ [/?e~ q /(e~ q ) - In (l - e^)] 
q 



(12) 



In Fig. @]we present the results of this entropy compari- 
son, showing how the initial temperature (T,) for the uniform 
system relates to the final temperature (Tf) in the deep lattice 
under the assumption of isentropic loading. We notice that 
Tf increases rapidly as Ti changes from zero. This arises be- 
cause for the gapless uniform system there is little energy cost 
to access low energy modes which contribute the finite initial 
entropy of the system. In contrast, in the deep lattice regime, 
a finite energy cost of order U must be supplied for the sys- 
tem to access excited states, so that k B T ~ U. Indeed, as is 
shown in Ref. (3, for k B Ti > 0.05S fl , and n — 1, the final 
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temperature scales linearly with U according to 

U ( E R \ 
T ' = ^( T * + °- m fcfj- (13) 

As U ~ V 3 / 4 , the temperature must increase with lattice 
depth, as is observed in Fig. 0] This is in contrast to the 
non-interacting results that show the temperature of the Bose 
system (in the low temperature regime where the results of 
this section hold) scales as eew ~ J, an d as noted earlier 
this parameter is exponentially suppressed with increasing lat- 
tice depth. Thus we find that the system has two competing 
behaviours: for low lattice depths, the system will be well- 
described by the non-interacting result and the temperature 
will decrease rapidly during loading, as observed in Figs. [3] 
(a) and (b). For deeper lattices, where tunneling between sites 
is small, the temperature is dominated by the energy gap of 
the excitation spectrum, i.e. U, which increases with increas- 
ing V. Including tunneling effects (i.e. finite J) somewhat 
suppresses the heating observed in Fig. [4] This can be quali- 
tatively understood in terms of the modifications that hopping 
makes to the eigenstate energies of the system. A non-zero 
value of J breaks the degeneracy of the energy gap U, lead- 
ing to a quasi-band whose width is proportional to J. As 
J increases (i.e when the lattice becomes shallower) the en- 
ergy of the lowest excited states decrease accordingly, while 
the ground state is only shifted by an amount proportional to 
J 2 /U. The lowest energy excitations then lie closer to the 
ground state and become accessible at lower temperatures. As 



a consequence, the entropy increases (and thus Tf decreases) 
with respect to the J — case. Of course in the deep lattice 
case, finite J corrections will become vanishingly small. 



V. CONCLUSION 

In this paper we have surveyed the physics of loading ultra- 
cold bosonic and fermionic atoms into optical lattices. Un- 
der the assumption that this loading is approximately adia- 
batic (isentropic) we have seen that there are regimes where 
the temperature of the system might be raised or lowered by 
the loading process. For bosons, the loading process can be 
used to reversibly condense the sample. For fermions, the 
Fermi energy sets a new energy scale, and for the case of a 
filled band no cooling regimes are available. We have exam- 
ined the effects of interactions on the Bose system, and seen 
that in the deep lattice limit the temperature of the system is 
proportional to the on-site interaction strength. 
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